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Optical properties of ensembles of three-level quantum emitters coupled to plas- 
monic systems are investigated employing a self-consistent model. It is shown 
that stimulated Raman adiabatic passage (STIRAP) technique can be success- 
fully adopted to control optical properties of hybrid materials with collective effects 
present and playing an important role in light-matter interactions. We consider a 
core-shell nanowire comprised of a silver core and a shell of coupled quantum emit- 
ters and utilize STIRAP scheme to control scattering efficiency of such a system 
in a frequency and spatial dependent manner. After the STIRAP induced popula- 
tion transfer to the final state takes place, the core-shell nanowire exhibits two sets 
of Rabi splittings with Fano lineshapes indicating strong interactions between two 
different atomic transitions driven by plasmon near- fields. 

PACS numbers: 32.80. Qk, 78.67.-n, 42.50. Ct 

INTRODUCTION 

The research field of plasmonics, while still expanding its applications in linear nano- 
optics [HE], is quickly advancing towards nonlinear phenomena [?H9]. Recently, it has been 
proposed to combine plasmonic systems with highly nonlinear media [10HT2j. Utilizing 
strongly inhomogeneous electromagnetic (EM) fields associated with the surface plasmon- 
polariton (SPP) resonance, one can achieve a significant spatial dependence of density of the 
conductive electrons in metals resulting in nonlinear phenomena such as second harmonic 
generation [131 E] • Extreme concentration of EM radiation was proposed to be utilized as a 
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catalyst to achieve lasing in nano-systems and has been recently implemented in experiments 
[To] . Other important applications of the light localization include surface enhanced Raman 
spectroscopy (SERS) [JB] and solar energy harvesting [TT| ITS]. 

A quickly growing field of hybrid materials is emerging [T91423] on the base of latest ad- 
vancements in nanoplasmonic science. Here one merges plasmonics with atomic and molec- 
ular physics considering systems comprised of quantum emitters and metal nano- structures. 
With current advances in chemistry and nano-manufacturing one is now able to couple molec- 
ular ensembles to plasmonic systems. Such systems exhibit a wide variety of new phenomena 
including new mixed molecular-plasmon states [21] and plasmon control of molecular energy 
redistribution [25] . 

It has long been realized that ideas of coherent control developed in quantum chemistry 
and physics [26] could be successfully applied to optically active nano-systems [27H3T] con- 
trolling electron transport [32], light pathways [33], and EM hot spots [34] . The ultimate 
goal in these investigations is to achieve control of optical properties of nano-structures. 

This paper explores ideas of stimulated Raman adiabatic passage (STIRAP) applied to 
ensembles of three-level atoms optically coupled to plasmonic systems. STIRAP is known 
to be based on the adiabatic population transfer within a single dressed state that does 
not include the dark transitional state thus minimizing spontaneous losses. The scheme 
has a variety of attractive modern applications from cooling internal degrees of freedom in 
molecules [33], to maximizing coherence between the initial and final states [SB], to manip- 
ulating dynamics in a multilevel system by making use of the Optimal Control Theory that 
reveals STIRAP type of control [37]. The goal of the paper is twofold - first, to analyze the 
efficiency of STIRAP technique in the ensemble of emitters where collective effects are taken 
into account in the framework of Maxwell-Liouville-von Neumann equations, and, second, 
to demonstrate the implementation of STIRAP as a tool to control scattering, reflection, 
and transmission properties of hybrid systems. As an example we consider a core-shell silver 
nanowire with resonantly coupled layer of three-level atoms. 

MODEL 

We consider electrodynamics of ensembles of quantum emitters in a self-consistent ap- 
proach. We solve the system of Maxwell's equations in time domain for electric, E, and 
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magnetic, H, fields. In spatial regions occupied by quantum emitters the Maxwell equations 
read 

BE „ -* dP 

"° aF = - v x E ' 

where £o and /xo and are the dielectric permittivity and the magnetic permeability of the 
free space, respectively, and P is the macroscopic polarization of a quantum medium. The 
latter is calculated using the mean-field approximation 

P = n a (d), (2) 

where (d) =Tr(pcf) is the expectation value of quantum emitter's dipole moment. The 
dynamics of the density matrix p satisfies the Liouville-von Neumann equation 

ihf t =[H,p\-iht{p)p, (3) 

here H is the Hamiltonian that describes the atom-EM field interaction. We assume that 
the relaxation processes are Markovian. The superoperator T accounts for the decay of the 
excited state and dephasing effects. 

In the mean-field approximation employed here it is assumed that the density matrix of 
the atomic ensemble is expressed as a product of density matrices of individual atoms (|3p 
driven by a local EM field (fil). In order to take into account dipole-dipole interactions of 
atoms within a single grid cell we follow Ref. and introduce Lorentz-Lorenz correction 
term for a local electric field according to 

— > 

E loc ^ = E + —, (4) 

6Eo 



where E is the solution of Maxwell's equations (Jlj) and macroscopic polarization is evaluated 
according to Eq. (pi). We performed several test simulations comparing results with and 
without local field correction term Eq. Q. It was found that frequency dependencies of 
observables, for example, the transmission coefficient is affected by Eq. Q at high atomic 
densities resulting in slight changes of resonant frequencies. If one is interested in qualitative 
analysis, such variations are not important. However, for the sake of completeness we use 
Eq. Q in all simulations below. 
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The system of equations ([T]) and ^ coupled via ^ with local field correction Q is 
solved on a parallel multiprocess cluster following the numerical algorithm discussed in Ref. 



RESULTS AND DISCUSSION 

First, we consider a A-system with the energy level diagram depicted in Fig. [l^, describing 
an atom in the ensemble of interacting emitters. An infinite in x and y dimensions and finite 
in z dimension quantum medium with the thickness of Az is driven by an incident linearly- 
polarized field. The vector of light polarization is along x axis and the light propagates 
along negative z axis as shown in Fig. [T]d. 

The Hamiltonian of a single three-level atom reads 



H = H - f2E } 



local 



Q Q -^local,^Ml2 ^ 

n fe, , E locai ^23 

huj 23 7T ^- 



(5) 



V 75 / 

where /zi 2 and /z 23 are transition dipoles (see Fig. JT^i). In all simulations we use /ii 2 = /i 2 3 = 2 

Debye. The Hamiltonian is written in the basis of states of angular momentum (|1), |3), 

|2)), the coefficients l/\/6 and l/-\/lO are the Clebsch-Gordan coefficients. 

As initial conditions we use an incident field in the form 

At (t — — — to) 2 
E x ; mc ,i = E 01 cos(u 12 (t ^ t )) exp( ^ ), (6) 



E X;iac o = E 02 cos(w 2 30 + to)) exp( 



Ar (t + f - t )- 
— - t )) exp( - 2 



'2 

where indexes 1 and 2 correspond to the pump field (to pump an atom from its ground state 
1 1) to the excited state |2)) and the Stokes field (to create a superposition of the |2) and the 
|3) states), with the latter preceding the former in time by Ar. Other parameters are the 
pulse duration r and the central time t . Their values are presented in the caption of Fig. 

m 

Our initial goal is to examine how STIRAP scheme for a single atom is affected by 
mutual EM interactions of atoms. In order to have reference data we first perform a series 
of simple calculations for a single atom case summarized in Fig. [2j The excitation pulse 
sequence is shown in the inset of Fig. [2^,. To determine optimal amplitudes of the Stokes 
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and pump pulses, E i, E 02 we calculate populations of atomic levels at the end of the pulse 
sequence scanning through their peak values E i and E 02 (here for simplicity we assume that 
E i = E 02 ). We chose the pulse duration T\ and r 2 to be such that 1/rj < (uj\ 2 — W23) = W13 
to resolve spectrally the splitting between the initial and final states. Fig. |2^i shows that the 
optimal STIRAP occurs at Eqi = Eq 2 = 1.1 x 10 9 V/m with the population of the target 
state |3) reaching 0.992 while the ground state population is 8 x 1CT 3 and the state |2) is 
nearly 0. In these calculations, the pulse duration r is chosen to be 1 ps, giving 1/r = 5a;i3. 
Note that longer pulse duration provides higher efficiency of population transfer up to 100%, 
our choice of the value of the r is dictated by computational reasons. Fig. [2]d shows time 
dynamics in a single atom under the optimal STIRAP condition confirming that the chosen 
scheme follows the conventional STIRAP. 

Next, we compare the results for a single atom with that obtained using one-dimensional 
self-consistent model for a layer of three-level atoms with thickness of Az = 200 nm. Panels 
(a) through (c) of Fig. [3] show the time dynamics of atomic populations averaged over 
the layer's volume at different atomic densities, n a , under the optimal STIRAP condition 
described above. At the density of n a = 1.5 x 10 27 m~ 3 the dynamics is nearly identical 
to that of a single atom. However at higher density of Ji a = 1.5 x 10 28 m~ 3 one can 
clearly see a noticeable difference. The final population of the target state degrades from 
its optimal value of 0.992 to 0.831 (Fig. [3^). It is interesting to note that even though 
STIRAP dynamics is significantly affected by strong coupling between atoms such that the 
intermediate state is populated during transitional times (Fig. [3|d), it is still negligibly small 
at the end of the excitation. One can examine a spatial distribution of both the ground and 
the target states at the end of STIRAP pulse sequence as shown in Fig. [3]d. Note that the 
incident field propagates from right to left (see Fig. [T]d for details). Atomic state populations 
exhibit spatial modulations with higher population transfer at the input side of the layer. 
This suggests that the pump and Stokes field carrier frequency gets modulated as the fields 
propagate through the medium. The degree of modulation is different for the Stokes and 
pump frequencies which results in a deviation from the two-photon resonance condition and, 
thus, in reduction of the efficiency of STIRAP population transfer. Another maximum of 
the target state population is seen near z = 37 nm with the corresponding minimum of the 
ground state population. STIRAP conditions, under which the sample is illuminated, are 
clearly affected by atom-atom interactions that alter incident field via spatially dependent 
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stimulated EM radiation. 

One may send two sets of STIRAP pulse sequences from both sides of the sample and 
examine the spatial dependence of atomic populations as shown in Fig. |4j Here we compare 
a single-sided excitation with a double-sided one. Several new features are noticed: a), the 
target state population becomes significantly higher in the center of the sample under sym- 
metric double-sided STIRAP condition; b). spatial modulations of atomic state populations 
are more pronounced. 

To illustrate one of many applications of the STIRAP scheme for optical control we 
perform a series of simulations calculating transmission, T, and reflection, R, coefficients of 
an atomic layer before and after STIRAP. Fig. [5] shows both sets of data. Before STIRAP 
pulse sequence partially inverts the atoms, the system has two reflection maxima near the 
transition frequency U12 with T exhibiting a wide minimum. The appearance of the second 
resonance in the reflection spectrum is a clear indication of a strong collective interaction of 
atoms in the layer that leads to appearance of new EM modes with frequencies other than 
Ui2- A number of these modes increases with the increase of atomic density [101 HT| . After 
the STIRAP population transfer to the final state, which is not coupled to the ground state, 
both R and T have extrema near U23 transition frequency. It should be noted that after 
STIRAP is applied, the transmission coefficient has an additional, small minimum at u^. 
This is due to the fact that not all atoms are inverted and their small fraction still produces 
that resonance. 

The second set of simulations is performed for a core-shell silver nanowire covered by 
a thin layer of three-level atoms as shown in Fig. [6^,. Here we perform two-dimensional 
simulations, assuming that the silver nanowire extends infinitely in z dimension. The system 
is excited by an incident electric field polarized along x axis that propagates along y axis. 

The energy level diagram of the system is created on an example of the alkali atom, e.g., 
Rb. Here, the electron transitions are considered that are induced between the hyperfine 
states of the 5 2 Si/2 and 5 2 Pi/2 electronic states within Dl line of 87 Rb (nuclear spin 1=3/2). 
For the 5 2 Si/2 state, F can take values 2 or 1, and for the Dl excited state 5 2 Pi/2, F is 
either 1 or 2. In our scheme, we choose the initial state to be F = 1 of the 5 2 5*i/2 state 
and the final state to be F = 2 of the 5 2 S'i/2 state. The lowest hyperfine sate F = 1 of the 
5 2 Pi/2 is the transitional excited state [12]. We consider the population initially to be in 
state F = 1 with the projection M = which we may prepare by optical pumping in the 
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presence of the constant magnetic field that removes the degeneracy of the magnetic sub- 
levels. Typical transition frequencies in plasmonic materials are on the order of 1 — 4 eV. In 
order to efficiently couple our model atom to SPP resonances we assume atomic transition 
frequencies to be on the same order as SPP modes. 

The Hamiltonian of a single three-level atom in the two-dimensional nanowire geometry 
shown in Fig. |6^i is 

/ 



H 




















































S1-M23 





2^3 


\/T0 




2^15 
C+^23 


2y/3 


2^15 





Q+[ll2 


S1-/112 


\ 


2^3 


2^3 


^-A«23 

Vio 







^+^23 


f^-M23 




2Vl5 


2Vl5 







f2+M23 




v^O 




hl0\2 










hoo 12 


J 



(7) 



V 

where Q± = Ei ocaljX ± iE\ ocg x ty . 

It is informative first to examine scattering efficiency of such a system in the linear 
regime, when |1) to |2) absorption line is dominant and is in resonance with a SPP mode 
of the silver nanowire. Under such conditions the localized SPP resonance in the scattering 
spectrum splits into two modes, upper and lower polaritons, the phenomenon known as Rabi 
splitting. The scattering spectrum for the bare silver nanowire and core-shell system are 
shown in Fig. [6Jd as functions of the incident frequency. The Rabi splitting is observed 
to increase with the density of the atomic ensemble. It is important to note that at high 
density of n a = 5 x 10 27 m -3 we observe additional resonance - the collective atom-plasmon 
mode. The latter has been detected several times in recent experiments [20J. The physics of 
this mode was recently scrutinized in Ref. [21], where it was demonstrated that this mode 
is due to plasmon induced dipole-dipole interactions between the quantum emitters. 

We apply the STIRAP excitation scheme with the parameters shown in the caption of 
Fig. [6] to the atomic layer covering silver nanowire as in Fig. [6^1. Time dynamics of atomic 
populations is shown in Fig. ^jp. We note that nearly perfect STIRAP, observed in a single 
atom case, is suppressed in the atomic layer. This is owing to several factors. First is a fast 
decoherence whose rate is chosen to be 10 13 s , which is one order of magnitude smaller 
than the peak Rabi frequency. (Here, we considered two main channels of decoherence, 
spontaneous emission and collisional dephasing.) The second factor reducing the efficiency 
of STIRAP relates to the spacial features of the sample: different locations of the atomic 
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layer are exposed to different EM fields due to spatially dependent strong local field enhance- 
ment near the surface of the silver core. Hence atoms in the shell are excited in a spatially 
dependent manner diminishing STIRAP. We performed additional simulations varying de- 
coherence rates (within femtosecond time scale) and the atomic density (on the order of 10 26 
to 10 28 m -3 ). In all simulations STIRAP was suppressed but still quite noticeable as in Fig. 
Eb. 

To utilize STIRAP technique as a possible control technique of optical properties of an 
atomic ensemble coupled to a plasmonic material, we calculate scattering intensity of the 
core-shell nanowire after STIRAP is complete. The results are shown in Fig. [6ji for the 
density of n a = 5 x 10 27 m~ 3 . Even though obviously not all atoms are inverted to the 
target state, the spectrum appreciably differs from the one corresponding to all atoms in 
the ground state. We note several important features: a), the Rabi splitting, that is 216 
meV before STIRAP, is reduced to 78 meV; b). the collective atom-plasmon mode is no 
longer seen; c). a new resonance near transitional atomic frequency U23 is observed due to 
the presence of inverted atoms. An additional resonant feature characteristic to inverted 
atoms is observed in the form of a tiny but distinct Rabi splitting in 3 meV. This is due to 
the fact that the SPP resonance of the silver nanowire is very broad and can be coupled to 
both |1) to |2) and |2) to |3) atomic transitions. Moreover the strong interaction between 
these transitions is also seen in the spectrum. Both resonances have evident Fano lineshapes, 
which indicates coherent interactions between the atoms that are in different states. 



CONCLUSIONS 



In the framework of Liouville von Neuman equation coupled to the Maxwell equations 
within the self-consistent approach, we demonstrated that STIRAP technique may be used to 
control optical properties of ensembles of quantum emitters coupled to plasmonic materials. 
The importance of such control parameter as the atomic density is emphasized. The results 
are obtained using self-consistent calculations where the STIRAP scheme is explored in one 
and two dimensions taking into account collective effects. It is shown that at low densities 
STIRAP scheme gives the result nearly identical to that for a single atom. Simulations 
at higher densities revealed the significance of collective interactions between atoms that 
eventually diminish the STIRAP control mechanism. When STIRAP scheme was applied to 
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hybrid nano-structures comprised of coupled three-level atoms and a silver nanowire, their 
scattering spectra manifest double Rabi splittings associated with two atomic transitions. 
It was also demonstrated that plasmon-polaritons induce strong interactions between these 
transitions leading to Fano lineshapes of scattering resonances. 

This work is partially supported by the National Science Foundation under Grant PHY- 
1205454 and grant from the Institute of Atomic, Molecular, and Optical Physics and the 
Smithsonian Center for Astrophysics. 
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FIG. 1: (Color online) One-dimensional STIRAP. 
atom with two dipole transitions indicated as blue 
simulations with a layer of atoms of the thickness 
along negative z direction and polarized vertically, 
set of parameters is used: W12 = 2.9 eV, W23 = 2.8 



(b) 




Panel (a): the energy diagram of a three-level 
and red arrows. Panel (b): a schematic setup of 
of Az exposed to the incident field propagating 
In all one-dimensional simulations the following 
eV, Az = 200 nm. 
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FIG. 2: (Color online) Single atom STIRAP. Panel (a): atomic state populations as functions of 
the incident field amplitude Eq (in V/m) for a single atom case. Ground state, |1), population 
is shown as a dashed black line, target state, |3), population is indicated as a dash-dotted red 
line, the solid blue line shows the intermediate state, |2), population. The inset depicts STIRAP 
scheme showing Stokes (solid black line) and pump (dashed red line) pulses as functions of time in 
ps. Panel (b): time dynamics of the atomic populations during STIRAP pulse sequence (here the 
color scheme is the same as in panel (a)). The parameters for STIRAP pulses (see eq. (|6|) are: 
E i = E 02 = E = 1.1 x 10 9 V/m, n = r 2 = 1 ps, At = 1.5 ps, t = 4 ps. 
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FIG. 3: (Color online) One-dimensional self-consistent calculations. Panel (a): spatially averaged 
ground state population as a function of time during STIRAP for a single atom case (solid black 
line), atomic layer at n a = 1.5 x 10 27 m -3 (dashed red line), and atomic layer at n a = 1.5 x 10 28 
m~ 3 (dash-dotted blue line). Panel (b): same as in panel (a) but for the intermediate state. Panel 
(c): same as in panels (a) and (b) but for the target state. Panel (d): spatial distribution of the 
ground state (solid black line) and target state (dashed red line) populations after STIRAP as 
functions of the coordinate z in nm at n a = 1.5 x 10 28 m . The pure dephasing rate is 0, the 
radiationless decay rate for both atomic transitions is 10 12 s . The transition frequencies for the 
atomic system are W21 = 2.8 eV, W32 = 2.9 eV. 



15 




z [nm] 

FIG. 4: One-dimensional self-consistent calculations. Spatial distribution of the ground state and 
target state populations after STIRAP as functions of the coordinate z in nm at n a = 1.5 x 10 28 
m~ 3 for two STIRAP excitation schemes: ground state population after a single-ended excitation 
- solid line, target state population after a single-ended excitation - dashed line, ground state 
population after the double-ended symmetric STIRAP - dotted line, target state population after 
the double-ended symmetric STIRAP - dash-dotted line. Other parameters are the same is in Fig. 
3, 
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FIG. 5: One-dimensional STIRAP control. Reflection (circles) and transmission (squares) coef- 
ficients as functions of incident frequency, co, in eV for the atomic layer at n a = 1.5 x 10 27 m~ 3 
before STIRAP (filled circles and squares) and after STIRAP (empty circles and squares). 
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time [ps] co [eV] 

FIG. 6: (Color online) Two-dimensional STIRAP control of core-shell nanowires. Panel (a): energy 
level diagram of a three-level atom in two dimensions. In all two-dimensional simulations = 3.61 
eV, W23 = 3.3 eV. The schematic setup of a core-shell nanowire is shown below the energy diagram. 
In simulations radii of the core and shell are set at R\ = 20 nm and i?2 = 35 nm, respectively. Panel 
(b): scattering intensity as a function of incident frequency, uj, in eV. The scattering spectrum of 
a bare silver wire is shown as a solid black line. Long-dashed red line shows data for core-shell 
wire at the atomic density of n a = 5 x 10 26 m~ 3 , short-dashed green line is for n a = 10 27 m~ 3 , 
and dash-dotted blue line is for n a = 5 x 10 m -3 . Panel (c): time dynamics of local atomic 
populations near surface of the sliver core during STIRAP. Ground state population - solid black 
line, intermediate state population - dash-dotted blue line, the target state population - dashed 
red line. Amplitudes of both Stokes and pump pulses is 3.4 x 10 9 V/m. The pure dephasing 
rate is 10 13 s _1 , the radiationless decay rate for both atomic transitions is 10 12 s _1 . Panel (d): 
same as in panel (b) but before STIRAP (dashed black line) and after STIRAP (solid red line) at 
n a = 5 x 10 27 m~ 3 . 



